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We present the analysis of the mean switching time and its standard deviation of 
short overdamped Josephson junctions, driven by a direct current and a periodic 
signal. The effect of noise enhanced stability is investigated. It is shown that 
fluctuations may both decrease and increase the switching time. 
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I. INTRODUCTION 

The investigation of nonlinear properties 
of Josephson junctions (JJs) is very impor- 
tant due to their broad apphcations in logic 
devices. Recently, a lot of attention was 
payed to Josephson logic devices with high 
damping. In papers [Rylyakov & Likharev, 
1999; Bunyk et al, 2001] a description 
and analysis of the entire system of single 
flux quantum logic elements are presented. 
The noise properties of systems consisting of 
many elements may be well understood from 
the noise properties of a single JJ. The pro- 
cesses occurring in such devices are based on 
a reproduction of quantum pulses due to 27r 
spasmodic change of the phase difference of 
the overdamped JJs. 

The use of high- Tc JJs creates many new 
problems, such as the thermally induced dig- 
ital errors and the effect of spreading of 
the switching speed due to thermal fluctu- 
ations. The analysis of influence of ther- 
mal fluctuations on different characteristics 
of JJs, such as the current-voltage charac- 
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teristic or the life time of superconductive 
state, has been made before on the basis of 
Langevin approach [Barone & Paterno, 1982; 
Likharev, 1986; Ortlepp & Uhlmann, 2004; 
Ortlepp & Uhlmann, 2005]. Recently exper- 
imental and theoretical work has been done 
on the mean switching time or mean escape 
time from the metastable state in periodically 
driven Josephson junctions [Yu & Han, 2003; 
Pankratov & Spagnolo, 2004; Peltonenet al, 
2006; Sun et al, 2007]. 

The process of switching of a single JJ 
with high damping under periodic driving 
was considered in [Pankratov & Spagnolo, 
2004]. In that paper the influence of the reso- 
nant activation (RA) and the noise enhanced 
stability (NES) phenomena in the accumula- 
tion or suppression of timing errors in rapid 
single flux quantum (RSFQ) devices was an- 
alyzed. Specifically an interval of frequen- 
cies was found where the switching time in- 
creases with increasing noise intensity. This 
effect is called noise enhanced stability [Man- 
tegna & Spagnolo, 1996; Agudov & Spagnolo, 
2001; Dubkov, Agudov & Spagnolo, 2004; Fi- 
asconaro et al, 2005; Spagnolo et al, 2007]. 
This paper presents the analysis of the ef- 
fects of thermal fluctuations on the switch- 
ing time for a single Josephson element with 



high damping, within the well-known resis- 
tive model of a Josephson junction [Likharev, 
1986]. The goal of this work is to investi- 
gate in detail how the NES phenomenon de- 
pends on the parameters of the periodical 
driving signal such as the amplitude and the 
frequency. 

II. GENERAL EQUATIONS AND 
STATEMENT OF THE PROBLEM 

The dynamics of a short over damped J J, 
under a current i{t) is given by the following 
Langevin equation 



.id(p{t) du{ip) 



Hit), 



dt dip 
u{ip) — 1 — cos{if) — i{t)(fi, 



(1) 



where (p is the phase difference of the or- 
der parameter [Barone & Paterno, 1982], 
u{(p) is the dimensionless potential profile 
(see Fig. 1.), Uc — 2eRNlc/^ is the charac- 
teristic frequency of the J J, Ic is the critical 
current, R]^ = Gjq \s the normal conduc- 
tivity of the JJ, e is the electron charge and 
h = h/2n, with h the Planck constant. Here 
i{t) = io + f{t) is the total current across 
the junction, zq is the constant bias current, 
f{t) is the driving signal, i{t) = I{t)/Ic, 
if it) = If{t)/Ic, and //(t) is the random com- 
ponent of the current. Because of thermal 
fluctuations, the random current may be rep- 
resented by white Gaussian noise 

{tf{t)) = 0,{tf{tUt + r)) = ^6{T), (2) 

where 7 = 2ekT/hIc = It/Ic is the dimen- 
sionless intensity of fluctuations, T is the 
temperature and k is the Boltzmann con- 
stant. 

Initially, the JJ is biased with a current 
across the junction smaller than the critical 
one, that is = {Iq/Ic) < 1, so as the initial 
condition we take the location of the phase in 
a potential minimum, (/?o = arcsin(io). A cur- 
rent signal f{t), such that i{t) > 1, switches 



therefore the junction into the resistive state. 
In Fig. 1 we show the periodical potential 
profile of the JJ and its extreme positions 
within which it varies in time. The switch- 
ing occurs not immediately, but at the later 
time, which is called the switching time. Due 
to the noise the switching time is a random 
quantity. We investigate therefore the mean 
switching time (MST) r and its standard de- 
viation (SD) a. As a driving signal we choose 
a sinusoidal signal fit) = Asm{u!t), where uj 
is the oscillation frequency and A is the signal 
amplitude. 

According to the definition [Malakhov & 
Pankratov, 2002] the mean switching time r 
and its standard deviation a are 

r = (i) = / tw{t)dt, {t"^) = / t^w{t)dt, 




(3) 

where w{t) = —dP{t)/dt, P{t) is the proba- 
bility to find if within the interval (— tt, tt). 
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FIG. 1 The potential profile of the JJ: the ex- 
treme positions during the periodical variations 
{% = —0.5 and i = 1.5), and the intermediate 
configuration (i = 0.5). The bias current and the 
signal amplitude are respectively: io = 0.5 and 
A = 1. The height of the barrier in the middle 
configuration of the potential profile is approxi- 
mately 0.7. 

In Fig. 2 the MST is shown versus the sig- 
nal frequency a; for = 0.5, /I = 1 and for 
different values of the noise intensity, namely 
7 = 0.05, 0.2, 0.5, 1. A frequency range, from 
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0.2 to 0.4, where the MST increases by in- 
creasing the noise intensity is clearly visible. 
In all the figures the MST r and the fre- 
quency cu are normalized to l/cUc and cUc, re- 
spectively. 
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FIG. 2 MST versus the signal frequency for 
different values of the noise intensity: 7 = 
0.05, 0.2, 0.5, 1. Here io = 0.5 and A = l. 

In Fig. 3 we report the same quantity 
of Fig. 2, that is the MST versus the sig- 
nal frequency, for a smaller signal ampli- 
tude A — 0.7 and different noise intensities, 
namely 7 = 0.005,0.01,0.05,0.2,0.5,1, and 
io = 0.5. Wc sec that, the frequency range 
where the NES effect is observed is reduced, 
and it is shifted towards smaller values of the 
driving frequency cu, namely from 0.1 to 0.2. 




0.03 0.1 0.3 1 



FIG. 3 MST versus the signal frequency for 
different values of the noise intensity: 7 = 
0.005,0.01,0.05,0.2,0.5,1. Here iq = 0.5 and 
A = 0.7. 



In both Figs. 2 and 3, the MST has a min- 
imum as a function of the driving frequency. 
This is the signature of the resonant acti- 
vation phenomenon, investigated in a JJ in 
Refs. [Yu & Han, 2003; Pankratov & Spag- 
nolo, 2004; Sun et aL, 2007]. 

III. RESULTS AND DISCUSSION 

The enhancement of the switching time 
may be considered in detail by plotting the 
dimensionless time MST function of 

the noise intensity for different vahics of the 
signal frequencies. Different behaviors of 
MST, depending on the values of the sig- 
nal frequency, occur for physical systems with 
metastable states [Agudov & Spagnolo, 2001; 
Dubkov, Agudov & Spagnolo, 2004]. Specif- 
ically a nonmonotonic behavior and a mono- 
tonic one may be observed. 

Let us start our consideration from small 
to large frequencies for the case io = 0.5 and 
A = 1, see Fig. 2. For frequencies smaller 
than 0.2 the character of the curves is sim- 
ilar: the MST decreases monotonically with 
increasing noise intensity. This case is pre- 
sented in Fig. 4 by two values of frequency: 
a; = 0.05 and 0.1. If the value of the noise in- 
tensity becomes greater or equal to the height 
of the barrier, which is approximately equal 
to 0.7 in the middle configuration (see Fig. 1), 
the particle does not see the fine structure 
of the potential. This is the reason why we 
restrict to 1 the values of the dimensionless 
noise intensity in the plots (Figs. 2-9). So, 
we can see that the MST decreases, with re- 
spect to the level of deterministic switching 
time (i.e. for 7 = 0), by increasing the noise 
intensity. The explanation of this behavior is 
the following. For small noise intensity and 
very low frequencies, the particle lies in the 
minimum for a long enough time due to the 
slow variation of the potential profile. The 
switching from the superconductive state to 
the resistive one does not occur until the po- 
tential barriers disappear almost completely. 
By increasing the noise intensity, the proba- 
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bility of the thermal activated switching in- 
creases and result the MST decreases. 

The next interval of frequencies which we 
consider is from 0.2 to, approximately, 0.52, 
because for such frequencies the behavior of 
the curves is similar. In this frequency range 
the NES effect appears, the MST (7) - curve 
behaviour becomes nonmonotonic in contrast 
with the previous case (a; < 0.2). In Fig. 4 
this behavior is shown for three frequency val- 
ues, namely u = 0.5, 0.45, and 0.4. We note 
that in these three curves, after the nonmono- 
tonic behavior, the MST increases again with 
the noise intensity 7 (for 7 > 1). Due to the 
increasing values of the signal frequency the 
escape process from the metastable state be- 
comes more rapid for low noise intensities and 
the MST for 7 — > is lower with respect to 
the previous frequency range {uj < 0.2). Due 
to the periodic variation of the potential pro- 
file and for noise intensities smaller than the 
barrier height, the particle starting from the 
minimum reaches a position near the top of 
the barrier. Then there is an optimum range 
of the system parameters, including the driv- 
ing parameters, for which the particle turns 
back at the metastable state because of the 
noise [Agudov & Spagnolo, 2001; Dubkov, 
Agudov & Spagnolo, 2004]. An enhancement 
of the lifetime of the metastable state pro- 
duces an enhancement of the mean switching 
time. By increasing the noise intensity the 
thermally activated escape increases too and 
the MST decreases. For a further increase of 
the noise, for intensity values greater than the 
barrier height (7 > 1), the potential profile 
"seen" by the particle is a linear one, which 
fluctuates between the extreme positions of 
Fig. 1. The particle can reach more easily po- 
sitions near the state at </? = and the MST, 
therefore, restarts to increase again with the 
noise intensity. 

The nonmonotonic behavior shown in 
Fig. 4 may be explained also by considering 
the competition of two factors in the above 
mentioned frequency range: the influence of 
the returning force from the lefthand side of 
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FIG. 4 MST versus noise intensity for seven 
different values of driving frequencies, namely 
Lo = 0.05,0.1,0.2,0.3,0.4,0.45,0.5. Here A = 1 
and iQ = 0.5. 



the potential profile and the increase of the 
variance of the phase position with increas- 
ing temperature. This explanation is related 
to the static case when the switching event 
starts because the value of the bias current 
is larger than the critical one [Malakhov & 
Pankratov, 1996]. By increasing the noise 
intensity the phase variance increases too. 
When the phase variance is not so large, the 
main part of the phase probability distribu- 
tion is located on the flat part near the point 
ip — arcsin(io), where the influence of the re- 
turning force is low. Thus, the influence of 
variance dominates, and the mean switching 
time increases with increasing of the noise in- 
tensity, because the particle stays more near 
the metastable state. Further increase of the 
phase variance, due to the noise, produces 
a larger probability distribution of the parti- 
cle, which reaches high slope positions of the 
potential profile. Then the returning force 
begins to affect from left to right because of 
the asymmetry induced by the potential pro- 
file and does not permit, in this frequency 
range, the probability distribution to move 
into the direction where the potential profile 
goes up quickly. Therefore, the probability 
distribution expands only into the right di- 
rection and the MST decreases because the 
particle escapes. Finally, the noise intensity 
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becomes so large that it is possible for parti- 
cle to move upstairs along the potential. The 
returning force is not able to hold the increas- 
ing of the variance anymore and the switching 
time increases again. 

By further increasing of the signal fre- 
quency {u > 0.5, for the parameter values 
used to obtain the curves in Fig. 4) a trap- 
ping phenomenon occurs. A threshold fre- 
quency uth exists which does not allow the 
particle to move to the next valley during 
one period of the signal. This means that 
for driving frequency u > Ufh, the particle 
is trapped within one period of the poten- 
tial profile and, as a consequence, the MST 
diverges (tends to an infinite value) with- 
out noise. The value of the threshold fre- 
quency increases with increasing the bias cur- 
rent and/or the maximal current across the 
junction [Agudov & Spagnolo, 2001; Dubkov, 
Agudov & Spagnolo, 2004]. The frequency 
dependence of the switching time, for zero 
noise intensity, is presented in the inset of 
Fig. 5. To be trapped the particle, starting at 
t = from the position ipQ (see Fig. 1) of the 
potential profile in the middle configuration 
{i = 0.5), should not reach the position at 
ip2, in the lowest configuration of the poten- 
tial profile {i = 1.5), within a quarter of the 
signal period. As a consequence the switch- 
ing event from the superconductive state to 
the resistive one doesn't occur. To estimate 
the threshold frequency we calculate, in ab- 
sence of noise, the time t spent by the particle 
to go from the point ipo to the point ip2, in the 
most favorable case, when the potential pro- 
file is fixed and is in the lowest configuration, 
with the current i = 1.5 (see Fig. 1). 

This time t, which is well approximated by 
one quarter of the threshold period Tt^, can 
be found from the solution of Eq. ([T]) without 
random current if and for a constant value 
of the total current i{t) exceeding the critical 
one [Malakhov & Pankratov, 1996] 

F((^2) - F{ipo) 



arctan 



i tan((y9/2) — 1 



•(4) 



From here the threshold frequency uoth = 
2'K/Tth — vr/2t. The values of ujth obtained 
from Eq. (j4]) are marked by the vertical lines 
in the inset of Fig. 5. 




FIG. 5 MST versus noise intensity for dif- 
ferent values of the driving frequency, namely 
uj = 0.55, 0.65, 0.7, 0.71, 0.72, which approach the 
threshold frequency uJth ~ 0.72 for a bias current 
io = 0.8 and a signal amplitude A = 0.7. In- 
set: the MST versus the driving frequency, from 
numerical simulations of Eq. ([1]) in the absence 
of noise (7 = 0), for three groups of values of A 
and zq. The different values oiujth, calculated by 
Eq- Qj are marked by vertical lines. 

Specifically we have: Uth = 0.48, for = 
0.5 and A = 1; ojth = 0.24, for = 0.5 and 
A = 0.7; uth = 0.69, for io = 0.8 and A = 0.7. 
As it is seen from the inset of Fig. 5, these 
values are very close to the exact values cal- 
culated by numerical simulations of Eq. ([1]) 
in the absence of noise (i/(t) = 0). The most 
prominent NES effect is observed when the 
frequency of the periodic signal is near the 
threshold frequency (see Fig. 5). For = 0.8 
and A = 0.7 the maximum MST is approxi- 
mately three times greater than its value at 

7 = io-\ 

In the next Fig. 6 the curves of MST are 
shown vs the noise intensity for five frequency 
values near and larger than the threshold Uth, 
that is for > 0.5 ~ Uth- Two distinct 
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FIG. 6 MST versus noise intensity for five 
values of the driving frequency, namely u = 
0.5, 0.6, 0.7, 0.8, 2. The bold curve gives the ex- 
act analytical behavior of MST for constant total 
current i = 0.5. Here A = 1, iq = 0.5. 

transient dynamical regimes are visible. For 
u = 0.5 we see a nonmonotonic behavior, 
with a finite value of r for 7^0. For fre- 
quencies larger than the threshold value Uth, 
all the curves are characterized by a mono- 
tonic divergent behavior in the limit of small 
noise intensity (7 0). In both dynam- 
ical regimes a minimum of MST is present 
for a noise intensity value of the order of the 
barrier height. By increasing the signal fre- 
quency the curves approach the behavior of 
the MST obtained with a fixed potential pro- 
file, corresponding to a constant total current 
i — 0.5. For high frequency values, in fact, 
the fiuctuations of the potential profile are so 
rapid that the potential "seen" by the parti- 
cle is the average potential, that is the mid- 
dle configuration of Fig. 1, with i = 0.5. In 
Fig. 6 the bold curve corresponds to the case 
of fixed potential profile with constant total 
current i = 0.5, and it was calculated by us- 
ing the exact analytical expression of MST, 
obtained in [Malakhov & Pankratov, 1996]. 

In the regime of frequencies above the 
threshold one the conditions of stochastic res- 
onance (SR) observation are fulfilled (see e.g. 
[Borromeo & Marchesoni, 2000]). SR is man- 
ifested when the MST does not exceed the pe- 
riod of signal and thus, the system response 




FIG. 7 MST and SD versus noise inten- 
sity for four driving frequencies, namely uj = 
0.1,0.3,0.45,0.5. The bold curve gives the 
asymptotic analytical behavior (7 << 1) of SD 
for constant total current i = 1.5. Here A = 1, 
io = 0.5. 

becomes periodic owing to noise (see Fig. 5, 
curve uj = 0.72 for which minimal MST is 
equal to ~ 6 and Fig. 6, curve u = 0.6, min- 
imal MST 10). 

In Fig. 7 the MST and its SD for four val- 
ues of the signal frequency are shown. We 
note that, rising from low to high frequencies 
near the threshold uJth, monotonic and non- 
monotonic behavior of the MST correspond 
to the same behavior of the SD, respectively. 
Moreover, in the limit of small noise inten- 
sity, all the SD behaviors show lower values 
with respect to the corresponding MST val- 
ues, but with different behavior depending 
on how much the frequency is close to the 
threshold Uth- For frequency u = 0.1, in fact, 
the SD increases slowly in such a way that 
it takes the same value of the MST for large 
noise intensity (7 = 1). But for greater fre- 
quency the noise intensity value, for which 
the two curves of MST and SD cross, de- 
creases. Specifically we have the following 
cross point values: 7 = 0.04 for u — 0.45 
and 7 = 0.002 for u = 0.5. This means 
that when the driving frequency is close to 
the threshold one the system becomes more 
unstable and even a small variation of the 
noise intensity has a strong influence on the 
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system. The curve, shown in Fig. 7 by the 
bold hne, is calculated for the case of the total 
current i = 1.5 and demonstrates a behavior 
of square root of the noise intensity (~ 1/7). 
This curve was obtained from the asymptotic 
analytical expansion of the SD in the small 
noise limit 7 << 1 and for a fixed poten- 
tial {u = 0) [Pankratov & Spagnolo, 2004; 
Gordeeva & Pankratov, 2006] . In the limit of 
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FIG. 8 SD versus noise intensity for five driv- 
ing frequencies, namely uj = 0.55,0.65,0.7,0.71, 
0.72. The bold straight line gives the asymptotic 
analytical behavior (7 << 1) of SD for constant 
total current i = 1.5. Here A = 0.7, iq = 0.8. 

— > the curves of the SD, obtained by our 
simulations, approach the asymptotic curve 
obtained for fixed potential. The behavior of 
the SD curve for = 0.1, obtained by nu- 
merical simulations of Eq. ([T]), is very close 
to the theoretical asymptotic curve {i = 1.5) 
not only in the low noise limit but also at 
larger values of 7 (7 = 0.1 — 1). This fea- 
ture of the SD(7) - curve behavior is shown 
in more detail in the next logarithmic plot 
(Fig. 8), which demonstrates the large rising 
of the SD for small noise intensities with in- 
creasing frequency. In Fig. 8 the curves of the 
SD are shown for the same frequency values 
and parameters used in Fig. 5 for the MST 
behaviors. Again we note that the cross point 
between the MST and SD curves, for frequen- 
cies close to the threshold one {u ~ 0.72 in 
this case) , is at very low noise intensities. For 
example for u = 0.71, the cross point is at 



7 = 0.15-10^'^. After the threshold frequency, 
we recover for the SD the same divergent be- 
havior of the MST, but with larger values. 
This fact confirms that the switching event in 
this region of parameters is due to the noise. 
It should be noted also that for large noise 
intensities (7 > 1) all graphs of the MST and 
SD collapse to the same curve, not depend- 
ing on the signal frequency (see Fig. 5 and 
Fig. 8). This means that the new regime of 
switching is unaffected by the signal and it is 
only due to the noise. 
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FIG. 9 MST versus noise intensity for two cases: 
(i) fixed bias current, namely io = 0.5 and A = 
1.05, 1.1, and (ii) fixed signal amplitude, namely 
^ = 1 and io = 0.5,0.55,0.6. Here lo = 0.5. 

In Fig. 9 the curves of MST vs 7 for differ- 
ent values of the bias current and signal am- 
plitude are shown. We consider two cases: (i) 
the bias current is fixed, namely io = 0.5 and 
A = 1.05, 1.1, and (ii) the signal amplitude 
is fixed, namely A = 1 and io = 0.5, 0.55, 0.6. 
We choose uj = 0.5 gnal frequency for 

which the NES effect is more pronounced. 
We see that changes of the bias current mod- 
ify the curves of MST more than changes of 
the signal amplitude. 

IV. CONCLUSIONS 

We analyzed the transient dynamics of a 
single overdamped Josephson junction driven 
by a periodic signal. The conditions of exis- 
tence of the NES effect are investigated. We 
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find that the enhancement of the switching 
time is larger for frequencies of the periodic 
signal close to the threshold frequency. In the 
region of the noise intensity values in which 
the mean switching time begins to increase, 
the standard deviation increases too but with 
different behavior depending on the proxim- 
ity of the signal frequency to the threshold 
Uth- The interval of frequencies, in which 
the NES effect is observed, depends on the 
maximal value of the total current across the 
junction and on the value of the bias cur- 
rent. There are two ways to change the value 
of the total current: by changing the am- 
plitude of the driving signal or the value of 
the bias current. We find that the switch- 
ing time changes more significantly due to 
the bias current variation, than due to the 
signal amplitude. We note that for small but 
nonzero values of the McCumber-Stewart pa- 
rameter the MST increases by several per- 
cent but qualitative character of the curves 
remains the same, as in the considered over- 
damped case. 
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